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Abstract. We present experimental features of identified particle production from 
nucleus-nucleus collisions at RHIC. These features reflect hadronization from a 
deconfined partonic matter whose particle formation scheme is distinctly different 
from fragmentation phenomenology in elementary collisions. Multi-parton dynamics, 
^ ' such as quark coalescences or recombinations, appear to be essential to explain the 

experimental measurements at the intermediate transverse momentum of 2-5 GeV/c. 
I Constituent quarks seem to be the dominant degrees of freedom at hadronization and 

gluon degrees of freedom are not explicitly manifested in the hadronization scheme. 
This physical scenario is consistent with recent Lattice QCD calculations that near the 



o 

^ . 

I critical temperature hadrons do not melt completely and quasi-particles may provide an 

■ effective degrees of freedom to describe the partonic matter. Heavy quark production 

^ , and other future experimental measurements to quantify the QCD properties of the 

produced matter at RHIC will be discussed. 

O 

-,— , , 1. Introduction 

X 

d • The advent of the Relativistic Heavy Ion Colhder (RHIC) at Brookhaven National 
Laboratory (BNL) has ushered in a new area in searches for the Quark Gluon Plasma 
(QGP). Recent measurements by all four RHIC experiments have established that a 
dense medium has been created in central Au-|-Au collisions at RHIC and the energy loss 
of partons leads to a suppression of high transverse momentum (p^) particles [H 12 El El • 
The QCD nature of the dense matter created at RHIC and whether the current 
experimental evidence proves the discovery of the QGP have been under debate within 
the heavy ion physics community IHl HI ■ We will attempt to address three specific 
questions: 1) features of hadronization and evidence for a color deconfined bulk partonic 
matter; 2) QCD properties of the matter at hadronization from experimental data 
and from Lattice QCD calculations; and 3) heavy quark production and other future 
experimental measurements which would further quantify the QCD properties at the 
phase boundary. 
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2. Features of Hadronization Prom Bulk Matter 

QCD color charges are confined and hadrons exist in color singlet state. The 
hadronization of quarks and gluons (partons) has not been understood based on QCD 
principles. Phenomenologically hadron formation in elementary e"*'e~ and nucleon- 
nucleon collisions has been described by fragmentation processes. In the pQCD domain 
physical processes are factorized into parton distribution functions, parton interaction 
processes and fragmentation functions for hadron formation. The fragmentation 
function is assumed to be universal and can be obtained from e^e~ collisions. A typical 
Feynmen-Field jB] fragmentation process involves a leading parton of momentum p, 
which fragments into a hadron of momentum ph whose properties are mostly determined 
by the leading parton. The fragmentation function is a function of variable z = "PhlVi 
where z is between and 1. Baryon production is found to be significantly suppressed 
compared to the production of light mesons: the baryon to pion ratio increases with z^ 
but never exceeds 20% [Oj. 

In the soft (low pj) particle production region, the pQCD framework and 
factorization break down - particles are believed not to be from fragmentation of partons. 
String fragmentation models were inspired by the QCD description of quark and anti- 
quark interaction. The Lund string model is one of the popular hadron formation models 
which has been successfully implemented in Monte Carlo description of e+e, p+p and 
nuclear collisions ^U]. In string fragmentation models the baryon production is also 
suppressed because the baryon formation requires the clustering of three quarks [TT]. 

Baryon production from nucleus-nucleus collisions, especially for multi-strange 
hyperons, has been measured to be much larger than theoretical model calculations. 
The hyperon production per number of participant pairs from nucleus-nucleus collisions 
at the SPS is significantly enhanced in comparison with the value from p+p collisions. 
The enhancement factor increases from A to f2 hyperons and with the collision centrality. 
The large production of strange hyperons cannot be described by string fragmentation 
models. In nucleus-nucleus collisions a large number of strings would overlap and lead 
to a much stronger string tension. A theoretical model to describe multiple strings as 
a rope has been proposed (12]. Strange hyperon production would be increased from 
ropes or strings with a larger tension. The large enhancement in the strange baryon 
production with respect to the number of participant scaling from p+p collisions has 
been proposed as evidence for the formation of the QGP at the SPS [13 . However, 
there is no direct experimental measurement linking the hyperon production with the 
QGP hadronization. Furthermore, theoretical calculations based on chiral properties of 
hadrons at high density fl] or strong color fields fH] can also provide an explanation 
for the large hyperon production. 

The increase in the baryon production from nucleus-nucleus collisions has become 
much more prominent at RHIC energies. Figure ^ shows the ratios of p/vr and KjKg 
from central Au+Au collisions at y^s^N = 130, 200 GeV measured by PHENIX |Tn] 
and STAR [HI dl CH] . The apparent difference between the STAR and PHENIX ratios 
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can be explained by the fact that the STAR A data include the electromagnetic decay 
contribution from TP and the PHENIX data are corrected for weak decay contributions. 
The STAR A data have been corrected for feeddown contributions from multi-strange 
hyperon decays. Some early data, for example [20] , were not included because these data 
were not corrected for weak decay contributions. The large baryon to meson ratio cannot 
be accommodated by the fragmentation scheme. The large ratio at the intermediate 
Pt region is the first indication that particle formation dynamics in nucleus-nucleus 
collisions at RHIC are distinctly different from the hadron formation mechanism via 
fragmentation processes in elementary e^e~ and nucleon-nucleon collisions. 
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Figure 1. Ratios of A to Ks from Au+Au and p+p collisions (STAR) and p to tt 
from Au+Au collisions (PHENIX) as a function of transverse momentum {pr)- The 
A data include contributions from E'^ decays and the p data are corrected for weak 
decays of hyperons. 



Measurements on nuclear modification factors and azimuthal angular anisotropy for 
identified particles have provided essential insight into hadronization in nucleus-nucleus 
collisions. The nuclear modification factor is defined as 

^ ^ [yield]^^ 

Nun X [yield]pv 

where Nun is the number of binary nucleon-nucleon collisions. The [yield]^"^ and 
[yieldY'^ are particle yields {d'^n/dprdy) from A-l-A and p-|-p collisions, respectively. The 
nuclear modification factor has also been defined using peripheral and central collisions 

^ ^ [yield/ NhinY"'''''''' 
[yield /NunY''''''^'''''"'^' 

A Raa or Rcp of unity implies that particle production from nucleus- nucleus collision is 
equivalent to a superposition of independent nucleon-nucleon colhsions. Hard scattering 
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processes within the pQCD framework are beheved to follow approximately binary 
scaling in the kinematic region where the nuclear shadowing of parton distribution 
function and the Cronin effect are not significant. Measurements of charged hadrons 
and neutral pions have revealed a strong suppression at high px region in central Au+Au 
collisions |2I1 1221 1231 121| ■ Recent d+ Au measurements [H |21 El El have demonstrated 
that the large suppression of high px particles in central Au+Au collisions is mainly 
due to energy loss, presumably of partons traversing the dense matter created in these 
collisions. 

Figure 121 shows the Rcp of K^, Ks, K*, 0, A, S and f2 as a function of pt from 
Au+Au collisions at ^/snn = 200 GeV measured by the STAR collaboration, where 
the Rcp ratio is derived from the most central 5% to the peripheral 40 — 60% collision 
centralities. The dashed line is the Rcp of charged hadrons for reference. In the low pt 
region soft particle production is dominated by the number of participant scaling. In 
the intermediate pr region of 2 to 5.5 GeV/c the pt dependence of Rcp falls into two 
groups, one for mesons and one for baryons. Despite the large mass differences between 
Ks and K*/(j), and between A and S little difference among the mesons and among 
the baryons has been observed within statistical errors. Particle dependence in the 
nuclear modification factor disappears only above a of 6 GeV/c, consistent with the 
expectation from conventional fragmentation processes. The unique meson and baryon 
dependence in the intermediate pt region indicates the onset of a production dynamics 
very different from both fragmentation at high p^ and hydrodynamic behavior at low 
Pt- 
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Figure 2. Rcp of K^, Ks, K*, (p, A, S and fl in comparison with charged hadron 
in dashed hne. Distinct meson and baryon groups are observed. 



The azimuthal angular distribution can be described by a Fourier expansion 
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where is the azimuthal angle of the particle, r is the reaction plane angle and V2 
has been called elliptic flow j2S]. Theoretical calculations indicate that V2 is generated 
mostly at the early stage of the nuclear collision. There are two dynamical mechanisms 
responsible for V2'- hydrodynamic and geometrical phase space of emitting source. In 
a non-central nucleus-nucleus collision the overlapping participants form an almond 
shaped particle emission source. The reaction plane is defined by the vectors x (impact 
parameter direction) and z (beam direction). The pressure gradient is greater and 
particles would experience larger hydrodynamic expansion along the short axis direction 
(in plane) than these along the long axis (out plane), resulting an ellipsoid in transverse 
momentum space in the final state. Particle production from a dense matter could lead 
to a surface emission pattern with higher particle density in the reaction plane than 
that out of the reaction plane. Figure El shows a schematic diagram for two scenarios of 
generating azimuthal angular anisotropy f 2 in non-central collisions. 
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Figure 3. Schematic diagram to sliow two dynamical origins of angular anisotropy, 
one from hydrodynamical expansion and the other from surface geometrical phase 
space. 



Figure m shows the V2 as a function of pt for vr, p, A and S from PHENIX [211 
and STAR [2Zj measurements. The angular anisotropy V2 reveals three salient features: 
1) particles exhibit hydrodynamic behavior in the low px region - a common expansion 
velocity may be established from the pressure of the system and the heavier the particle 
the larger the momentum from the hydrodynamic motion leading to a decreasing 
ordering of f 2 from vr, to and p for a given pt] 2) f 2 values do not depend strongly on 
Pt at the intermediate pt region in contrast to the strong p^ dependence in the yield of 
particles; 3) the saturated f 2 values for baryons are higher than those for mesons and 
there is a distinct grouping among mesons and baryons. 

The absence of a strong pt dependence at the intermediate pt region is an intriguing 
phenomenon for the angular anisotropy of particle emission. A very plausible scenario 
would be to relate the V2 in this pt region with the geometrical shape of the emitting 
particle source. A surface emission pattern from the almond shaped participant volume 
would naturally lead to a saturation of f2. A surface emission scenario is possible if 
particles are produced from a dense bulk matter and the hadronization duration is 
relatively short. 
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Figure 4. Azimuthal angular anisotropy V2 as a function oipT for identified particles. 
The hydrodynamic calculation is provided by P. Huovinen et al |28|. 

Parton energy loss in the dense medium created in nucleus-nucleus collisions was 
proposed as a possible mechanism for generating an angular anisotropy V2- High energy 
partons are quenched inside the dense medium and lead to an effective particle emission 
from a shell area of the participating volume [5^. Gyulassy and Vitev proposed the 
parton energy scenario as a mechanism to explain the deviation of f 2 from hydrodynamic 
behavior at moderate pt ISD] • While this scenario may be important for the considerable 
V2 magnitude for charged hadrons at a px greater than 6 GeV/c or so measured by 
STAR [HI], it cannot explain the particle dependence in both the nuclear modification 
factor and the angular anisotropy f 2 at the intermediate px region. Within the parton 
energy loss scenario the larger f 2 of baryons implies a higher energy loss than that of 
mesons; the larger nuclear modification factor of baryons, however, is only consistent 
with a smaller energy loss than that of mesons. 

Shuryak [32] has pointed out that the magnitude of the measured v 2 is significantly 
larger than what can be accommodated based on particle emission from a geometrical 
ellipsoid source within an energy loss scenario. Using a more realistic Wood-Saxon 
description of the colliding nuclei for the ellipsoid source the predicted theoretical V2 
is much smaller than that from a hard-sphere model of the colliding nuclei, leading to 
a greater discrepancy between the measurement and the theoretical expectation. The 
magnitude and the particle dependence of f 2 at the intermediate px region cannot have 
a dynamical origin either from hydrodynamic flow or from parton energy loss alone. 

An empirical constituent quark number (n) scaling has also been observed |33j . 
Figure El presents the V2/n as a function oipx/n for tt, K, p, A and S from Au+Au 
collisions, where the line is a fit to the data points excluding the vr data. The bottom 
panel shows the ratio of data points to the fit. At the intermediate px region ( 0.6 < 
Px/n < 2 GeV/c) all meson and baryon data points fall onto an uniform curve. In the 
lower Px region there is a particle dependent deviation from the fit curve presumably due 
to the hydrodynamic flow as a function of mass of particles. The vr data are significantly 
above the V2 of other mesons. The large fraction of resonance decay contribution to the 
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TT yield is known to enhance the f 2 of tt at a given pt IHH ISH] • 
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Figure 5. Azimuthal angular anisotropy V2/n as a function of px/n for identified 
particles where n is the number of constituent quarks. The line is a fit to the data 
points excluding the tt data and the bottom panel shows the deviation from the fit 
line. 

The quark number scahng hints at dynamics for the hadron formation at the 
intermediate px which is very different from fragmentation at high px and from thermal 
statistical descriptions at the low pt region. The scaled azimuthal angular anisotropy 
(f2/n) may be interpreted as the constituent quark anisotropy just prior to the hadron 
formation. The essence of this hadron formation dynamics is very different from the 
fragmentation picture where the leading parton plays a dominant role in determining 
properties of hadron. The measured features for intermediate pt hadrons produced 
at RHIC require that all quark ingredients (n=2 for mesons and 3 for baryons) play 
an approximately equal role in hadron formation and that the hadron properties are 
determined by the sum of partons. 

Recent theoretical models of quark coalescence [SEl IHII and recombination jHEl tHH] 
all have the essence of multi-parton dynamics for hadron formation despite significant 
differences in details. The recombination models provided a satisfactory description of 
the particle yields, in particular, the large production of baryons in the intermediate pr 
region. The formation of a dense partonic system provides a parton density dependent 
increase in baryon yield as a function of collision centrality through the coalescence 
mechanism. More importantly the particle dependence in f 2 requires a V2 distribution 
at the constituent quark level. The establishment of the quark collectivity V2 would be 
an unambiguous signal for deconfinement in bulk partonic matter. 
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3. Color Deconfined Matter and its QCD Properties at Hadronization 

The fact that for hadron formation at intermediate pt all quark constituents must 
contribute almost equally to hadron properties does not depend on details of theoretical 
models. Recombination |HHl inH] or coalescence models jHEl IHZ| provide a useful 
theoretical framework to derive quantitative quark level properties. However, these 
models do not constitute unambiguous evidence for deconfined partonic matter. The 
recombination model developed by Hwa and Yang [10], for example, works for p+p 
and d+Au collisions as well. These models provide an alternative hadronization scheme 
which, though different from fragmentation, has been demonstrated to be more suitable 
for hadron formation at the intermediate px region from nucleus-nucleus collisions at 
RHIC. The unique evidence for a color deconfined quark matter comes from the fact 
that constituent quarks must have a collective V2 distribution in order to explain the 
measured particle dependence of azimuthal angular anisotropy. 

The particle dependences in the nuclear modification factor and in the angular 
anisotropy V2 have demonstrated that two distinct groups: mesons and baryons, 
characterize the particle dependence at the intermediate pt region. The quark number 
scaling in the V2 indicates that constituent quarks have developed a collective angular 
anisotropy distribution f 2 in transverse momentum space by the time of hadronization 
and hadrons are formed through coalescence or recombination of constituent quarks. 
The dominant degrees of freedom at hadronization seem to be in the constituent quarks. 
Based on nuclear parton distribution functions and quark versus gluon interaction cross 
sections, the initial conditions of nucleus-nucleus collisions at mid-rapidity at RHIC 
are expected to be dominated by gluon degrees of freedom. We do not know how 
the colliding system evolves from the initially gluon-dominated matter to a constituent 
quark dominated system just prior to hadronization. The dominance of the constituent 
quark degrees of freedom and the collective f 2 distribution at the constituent quark level 
are the most direct experimental evidence for a color deconfined matter, despite the lack 
of our understanding on the dynamics of nuclear matter evolution. The saturation of 
V2 can be related to the geometrical shape of the particle emitting source at the time 
of hadron formation, which supports the notion of bulk partonic matter of constituent 
quark degrees of freedom. 

The emerging physical picture implies that constituent quarks are the dominant 
degrees of freedom at the boundary of quark-hadron transition. The gluon degrees 
of freedom are not explicitly manifested in the hadron formation or in characterizing 
the hadron properties. It appears that the nuclear matter created in nucleus-nucleus 
collisions at RHIC would evolve naturally from a partonic gluon-dominated initial state 
towards a constituent quark or quasi-hadron state at the time of hadron formation. 
These experimental observations may provide useful insights on the QCD properties of 
the quark matter near T^.. Recent Lattice QCD (LQCD) calculations indicate that 
spectral functions of pseudoscalar and vector mesons have non-trivial shapes at a 
temperature above the critical T^. In particular, heavy quarkonia such as J/?/' may 
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survive at a temperature above 1.6Tc |42[ l43j. Other theoretical calculations, for example 
reference jH], have also invoked the notion of quasi-hadrons to describe properties of the 
dense matter created at RHIC. It is a critical step to firmly connect the experimental 
insights on the properties of the quark matter at the boundary of hadronization with 
the LQCD calculations. The disappearance of the gluon degrees of freedom from the 
initial state and the emergence of constituent quarks at hadronization are some of the 
critical conceptual questions to be addressed. These questions will have a direct bearing 
on the nature of the quark hadron phase transition for which there is as yet no direct 
experimental evidence. 

4. Charm Production Cross Section at RHIC 

Heavy quarks are produced mostly in initial parton scattering or during the very high 
temperature phase of the collision. Therefore, heavy quark measurement can probe the 
initial parton flux, the dynamical evolution and the quark energy loss in dense medium. 
If heavy quarks are found to participate in the collective motion of the medium (radial 
or elliptic flow), this will lend further confirmation for parton collectivity. Features of 
charm meson production will also bear signatures of the partonic matter at the phase 
boundary. Both the STAR and PHENIX collaborations at RHIC have been pursuing 
vigorous heavy quark physics programs both in analyses of current data and in future 
detector upgrade plans to provide better capabilities for heavy quark measurements. 

The total charm quark pair production cross section (cTcc) is an important constraint 
on the collision dynamics and the heavy quark evolution. Both STAR and PHENIX have 
presented results on the cXcc measurement of p+p collisions from charm semi-leptonic 
decays. In addition, STAR has also derived an equivalent o"cc for p+p collisions based on 
direct reconstruction of hadronic decays of charm mesons from d+Au collisions. FigurelHl 
shows the px spectra of electron and from STAR. The PHENIX preliminary non- 
photonic electron data are represented by the fitted line with a reported cXcc measurement 
of 709 ± 85(stat)+332 — 281 (sys) [ih |33|. STAR has measured a cross section of 
1300 ±200 ±400 yu6 and a mean transverse momentum for of 1.32 ±0.08 GeV/c from 
direct reconstruction [IH]- A next-to-leading order pQCD calculation of the charm 
quark production cross section jlTj has yielded 300 to 450 /xfe, significantly below the 
STAR measurement and at the lower end of the PHENIX range of uncertainty. 

Several comments on the cross section measurements are in order. When measuring 
the charm cross section through semi-leptonic decay, the quality of the electron data 
below apT of 1 GeV/c or so is significantly deteriorated because of a large combinatorial 
background. An electron of ~ 1 GeV/c typically comes from the decay of a D 
meson with p-r ~ 2 GeV/c, which is significantly beyond the average pt of 1.32 ± 0.08 
GeV/c reported by STAR. Therefore, one has to extrapolate to the low pt region by 
over a factor of two to obtain the total charm cross section. Such an extrapolation is 
often model dependent and has a large uncertainty. The semi-leptonic decay branching 
ratios for D^, D*, and Ds are different. The electron yield from decays of these 
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Figure 6. pT distribution of D' 
decays of charm mesons. 



mesons and non-photonic electrons from semi-leptonic 



D mesons depends on the relative yield which is another important contribution to the 
uncertainties of the charm production cross section derived from electron measurement. 
The direct reconstruction of the D decay kinematics provides a broad coverage of pt and 
does not suffer from the same uncertainties as the leptonic decay electron measurement. 
However, present STAR measurements of D meson yields using event-mixing methods 
from TPC tracks suffer from limited statistics. A future vertex detector upgrade capable 
of measuring the D decay vertex displacement is essential for both STAR and PHENIX 
heavy flavor physics programs. 

Figure [7| presents the STAR preliminary transverse momentum spectrum of D 
mesons from d-|-Au collisions normalized by the number of binary collisions, where 
the Pt shapes of D* and are assumed to be the same as that of jlSl- The shape 
of the Pt distribution coincides with the bare charm quark pt distribution from the 
Fixed-Order- Next-Leading-Log (FONLL) pQCD calculation from M. Cacciari et al |49j . 
If a fragmentation function such as the Lund fragmentation scheme jSO] or the Peterson 
function [SI] is introduced for D meson production, the resulting pt distribution will 
be significantly below the measurement at the high pt region. This observation raises 
an outstanding question regarding the Pt distribution and the formation mechanism of 
D mesons in hadro-production. Recently a /ct factorization scheme has been found to 
significantly change the D meson pt distribution from nuclear collisions as well 1^21 • 

The fact that the D meson pt distribution can be better described by the pt of 
bare charm quarks from pQCD calculations has been observed in previous fixed target 
experiments jSHj- With a fragmentation function such as the Peterson function the D 
meson pt distribution is too steep to explain the measured pt distribution. In order 
to match the calculation with the data one has to introduce a kT kick to the parton 
distribution. The scale of < k^ > is on the order of 1 (GeV/c)^, much larger than the 
typical Aqcd scale for strong interaction. Furthermore, the Feynman xp of the D meson 
distribution was also found to coincide with the bare charm quark xp distribution |53j . 



Salient Features of Hadronization, Deconfinement and Heavy Quark Probes at RHICll 




Figure 7. STAR preliminary measurement of pt distribution of D mesons from 
d+Au collisions normalized by the number of binary collisions. The shape of the pt 
distribution is compared to pQCD FONLL calculations. 



The fragmentation function will have a large impact on the xp distribution from bare 
charm quark to D meson, which cannot be negated by introducing any longitudinal 
boost of reasonable scale as in the case for kx kick in the transverse momentum direction. 

The transverse momentum distribution of particles produced at RHIC energies 
are considerably flatter than those at lower energies. The k^ kick scheme does not 
change the shape of the pt distribution significantly. The STAR measurement of the D 
meson pj- distribution suggests either that the charm quark fragmentation may be close 
to a delta function or that a charm formation mechanism such as the recombination 
model may be important in hadro-production. The recombination model takes a charm 
quark and combines it with a light quark, presumably of low pt, to form a charmed 
meson. Therefore, the pt of the meson is not significantly different from that of the 
bare charm quark. Measurements of charmed meson production provide unique probes 
for the hadron formation dynamics and for the transport dynamics of heavy quarks in 
the dense nuclear medium produced at RHIC energies. 

5. Recombination Mechanism and Heavy Quark Flovi^ and Energy Loss 

R.J. Fries et al j38j has argued that recombination dominates over the fragmentation 
production mechanism when the parton spectrum is a thermal distribution. The 
resulting hadron momentum distribution exhibits the same thermal distribution as the 
constituents. Only at sufficiently large pr, where the underlying parton distribution is 
a power-law, does fragmentation overtake recombination. The recombination process 
is particularly effective for baryon production which recombines three readily available 
quarks while in fragmentation process the three quark production is suppressed. Since 
the low Pt distributions of hadrons from p+p, d+A and A+A collisions at high energies 
such as RHIC are not drastically different quanlitatively, the argument by R.J. Fries 
et al would imply that the recombination process should play a unique role at the 
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intermediate pt region for all these collisions if we interpret the soft particle production 
as from an underlying thermal parton distribution. Recombination model calculations 
by Hwa and Yang can fit Au+Au and d+Au data from RHIC reasonably well |4()j . 
The fact that phenomenologically an underlying parton distribution is used for the 
recombination model does not necessarily mean the formation of partonic matter. 

The experimental evidence of bulk partonic matter formation and its collectivity 
comes from the elliptic flow V2 of partons based on constituent quark scaling in the V2 
measurement of mesons and baryons. Therefore, one critical test of the idea is to measure 
many more particles and resonances over a much broader region. In particular, we 
need to push the pt coverage above 6 GeV/c or so to reach the fragmentation region. 
The meson and baryon difference in the nuclear modification factor and the angular 
anisotropy V2 is expected to vanish in the fragmentation region. 

STAR preliminary measurement of S and its comparison with A indicated that the 
strange quark dynamics in nuclear modification factor and the V2 at the intermediate 
Pt region is similar to these of light quarks: the constituent strange quarks have a 
similar V2 as light up and down quarks [HI]. The measurement of nuclear modification 
factor and f 2 of meson and Q baryon covering the full Pt range from recombination to 
fragmentation region in nucleus-nucleus collisions at RHIC will further test the strange 
quark dynamics in the hadronization of bulk partonic matter. Possible variation may 
come from the particle formation differences because the meson is a vector meson and 
fl belongs to decuplet while what have been measured are for pseudo scalar mesons and 
octet baryons. Resonance measurement has also been proposed to test the production 
dynamics ^E^. 

Charm quark transport dynamics in dense nuclear medium will provide unique 
probes to the QCD properties of the medium. If the initial temperature is very high, 
T ~ 500 MeV or higher, the yield of total charm quarks can also be increased through 
thermal gluon-gluon scatterings. Possible suppression of charm mesons at high pt 
will test the energy loss dynamics of charm quark propagation in a QCD medium. 
Theoretical calculations have predicted a reduced medium induced energy loss for heavy 
quarks and the high p^ suppression of charmed hadrons should not be as strong as light 
hadrons fS^^ EZI- Observation of heavy quark hydrodynamic flow would indicate that 
heavy quarks, once created in the initial state, must have participated in the partonic 
hydrodynamic evolution over a sufficiently long period of time to reach a substantial flow 
magnitude. This would be a unique probe for the early stage of a partonic phase jHSj- 
Figure |H1 shows preliminary STAR jSH] and PHENIX [^01 measurements of the f 2 for 
electrons from charm leptonic decays, which have been demonstrated to be closely 
correlated with charmed meson V2 |35j . 

6. Summary 

Characteristics in hadrons below a of 5 GeV/c or so from nucleus- nucleus collisions 
at RHIC have shown distinct features that are drastically different from fragmentation 
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Figure 8. Preliminary measurement of non-photonic electron V2 as a function of pt 
from PHENIX and STAR. The curves are from a calculation by V. Greco et al |58) . 



processes in elementary collisions. A salient feature of meson and baryon dependence 
has been observed in the nuclear modification factor and the angular anisotropy f 2 of tt, 
K"^, Ks, K*, p, A and S particles at the intermediate pt of 2-5 GeV/c. A constituent 
quark number scaling has been observed for the V2 measurement. These experimental 
measurements suggest that at the hadronization moment the effective degrees of freedom 
are the constituent quarks; the constituent quarks have developed a collective V2 as a 
function of pt] and the hadron formation at the intermediate pt is likely through multi- 
parton dynamics such as recombination or coalescence process. 

The physical picture emerging from the experimental measurements complements 
the Lattice QCD results. Spectral function calculations have indicated that hadrons, 
particularly heavy quarkonia, do not melt completely at critical temperature. It 
appears plausible that the constituent quark degrees of freedom or quasi-particles play a 
dominant role at the hadronization of bulk partonic matter though further confirmation 
of the picture from LQCD is needed. Despite intriguing experimental observations of 
hadronization from a deconfined bulk partonic matter, signatures for a quark-hadron 
phase transition remain elusive. 

Heavy quark production and its transport dynamics in dense nuclear medium 
probe QCD properties of the matter. The charm quark flow measurement will provide 
a significant insight on recombination or coalescence hadronization mechanism and 
partonic collectivity of the dense matter. Future detector upgrades from STAR and 
PHENIX will greatly enhance their heavy quark measurement capabilities at RHIC. 
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